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ABSTRACT 

A computed list of H"2 16 infra-red transition frequencies and intensities is presented. 
The list, BT2, was produced using a discrete variable representation two-step approach 
for solving the rotation- vibration nuclear motions. It is the most complete water line 
list in existence, comprising over 500 million transitions (65% more than any other 
list) and it is also the most accurate (over 90% of all known experimental energy levels 
are within 0.3 cm -1 of the BT2 values). Its accuracy has been confirmed by extensive 
testing against astronomical and laboratory data. 

The line list has been used to identify individual water lines in a variety of ob- 
jects including: comets, sunspots, a brown dwarf and the nova-like object V838 Mon. 
Comparison of the observed intensities with those generated by BT2 enables physical 
values to be derived for these objects. The line list can also be used to provide an opac- 
ity for models of the atmospheres of M-dwarf stars and assign previously unknown 
water lines in laboratory spectra. 
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1 INTRODUCTION 

Water is the most abundant molecule in the universe after 
H2 and CO. It is present in many astrophysical environments 
including the atmospheres of: M dwarfs (Allard et al.1994), 
brown dwarfs (Allard et al. 1996), K and M giants and su- 
pergiants (Jennings and Sada 1998; Ryde et al. 2002; Tsuji 
2001) and oxygen-rich AGB stars (Barlow et al. 1996). It 
occurs in: sun-spots (Wallace et al. 1992; Polyansky et al. 
1997), nova outflows (Banerjee et al. 2005), Mira variables 
(Hinkle and Barnes 1979), T Tauri eruptive variables (Shiba 
et al. 1993), dark molecular clouds (Gensheimer et al. 1996), 
young stellar objects (Carr et al. 2004), comets (Mumma 
et al. 1996; Dello Russo et al. 2000), the ISM (Cernicharo 

1994) , masers (Cheung et al. 1969; Gonzalez- Alphonso et al. 

1995) and planetary atmospheres. An accurate water line 
list is thus essential for interpreting spectra from all of these 
sources and in modelling stellar atmospheres at tempera- 
tures up to 4,000 K. 

The importance of water has given rise to many labo- 
ratory investigations of its spectrum. Ludwig 1971; Camy- 
Peyret et al. 1977; Bernath 1996 all investigated hot water 
line positions. However, techincal problems and the huge 
number of transitions (many of which appear blended) mean 
that only in the region of 80,000 (out of a total of more than 
a billion) transitions are known experimentally and there 
are few hot water lines for which intensities have been de- 
termined. 

The spectrum of water, which extends over a wide wave- 
length range from millimetre to near ultra-violet, is due to 



quantised changes in the rotation-vibration energy of the 
atomic nuclei moving in the electronic potential well. Essen- 
tially, the water molecule is only able to absorb or emit in 
its ground electronic state as the energy of the first stable 
excited electronic state is above the dissociation energy. In 
practice, some emissions do occur from short-lived excited 
electronic states. H2O electronic transitions from diffuse in- 
terstellar clouds provide an example (Smith et al. 1981), 
but these transitions, which occur in the ultra-violet, can be 
disregarded in almost all other situations. 

Water is a triatomic asymmetric top molecule. Its 
rotation- vibration spectrum is more complicated than those 
of most other triatomic molecules. In common with all non- 
linear triatomics, H2O has six degrees of internal freedom 
(three of rotation and three of vibration). However, the light- 
ness of the hydrogen atoms means that the rotation con- 
stants are large, and this gives rise to an open spectrum 
that extends over a wide frequency range. Moreover, the 
'floppy' nature of the molecule means that the movement 
of the hydrogen atoms is generally anharmonic and con- 
sequently transitions involving changes of more than one 
vibrational quanta often occur. Also, since many of the vi- 
brational frequencies are nearly resonant with other frequen- 
cies, it is common for vibrational bands to overlap and for 
states to interact in ways that cannot easily be predicted 
by perturbation theory but are amenable to a variational 
approach. 

The importance of the H2O molecule in astronomy and 
the complexity of its spectrum have created a great deal of 
interest in the possibility of generating the spectrum syn- 
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thetically. Previous synthetic line lists include: MT (Miller 
et al. 1994), VT1 (Viti et al. 1997) VT2 (Viti Ph.D.thesis 
1997), PS, otherwise called AMES, (Partridge & Schwenke 
1997) and SCAN (j0rgensen et al. 2001). All previous syn- 
thetic line lists have suffered from a number of problems that 
are discussed below. The most successful attempts have em- 
ployed similar variational nuclear motion procedures to that 
used in producing the BT2 water line list. However, for rea- 
sons detailed in the next section, none of the earlier lists is 
considered to be satisfactory. We have addressed these prob- 
lems and consequently the BT2 line list is an accurate tool 
for astronomers working in a variety of fields. 



Table 1. Symmetry Blocks 
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q is the vibrational basis symmetry and p is the rotational par- 
ity. These are labelled as either symmetric (e), or antisymmetric 
(o) states. O/P = Ortho/Para. The 'Code' is the notation for 
symmetry used in the Levels File (Table 2). 



2 BACKGROUND TO THE CALCULATIONS 
2.1 Variational Techniques 

Variational techniques represent the best approach to solv- 
ing the the nuclear motion problem (Tennyson 1992), and 
they are examined in detail in Bacic and Light (1989). Here 
we use a discrete variable representation (DVR). 

In a DVR, the wavefunctions are defined by a complete 
set of weighted, orthogonal grid points, each wavefunction 
having a different set of weightings. This method is capable 
of generating accurate solutions, the accuracy being deter- 
mined by the the number and appropriateness of the points. 
It is efficient for a large number of situations. The DVR 
approach has the advantage that the potential matrix el- 
ements are diagonal, and hence are easily evaluated. Even 
more important is the fact that the dipoles can be reduced 
to a similar form (Tennyson et al. 2004). 

The current work is not unique in employing a DVR 
approach to solve the nuclear motion problem for water (see 
for example, Viti et al. 1997). However, improved physics, 
in the form of a highly accurate potential energy surface 
(PES), the methodology embodied in the DVR3D program 
suite and increases in computational power have made it 
possible to produce a line list that is more complete and 
more accurate than any previous list, even those employing 
similar methodolodies. 

The DVR3D approach generates four rotation-vibration 
symmetry blocks for the H2O molecule, which we label ee, 
eo, oe and 00. The first e/o term is the vibrational basis 
symmetry 'q' and the second, e/o is the standard quantum 
number 'p', the rotational parity. These symmetry blocks 
are not the same as the C2«(M) symmetry blocks r r „: Ai, 
A2, Bi, B2, but are related to them (see standard texts on 
molecular spectroscopy, such as Bunker & Jensen 2005). 

The nuclear permutation operation in which the two 
identical protons comprising the hydrogen nuclei are inter- 
changed gives rise two two dicrete states of the molecule: or- 
tho (O) and para (P).The nuclear spins may couple symmet- 
rically or antisymmetrically. The antisymmetric coupling 
that gives rise to the O form of the H2O molecule is triply 
degenerate, whilst the symmetric coupling, that gives rise to 
the P form of the molecule is undegenerate. 

For any J, apart from J=0, there are two ortho and 
two para symmetry blocks (in the case of J=0 there is one 
O and one P block) There are two possible arrangements, 
depending on whether J is odd or even and these are detailed 
in Table 1. 



The difference in degeneracies of the O and P states 
impacts on the partition function of the molecule as well as 
on line intensity. Also, the fact that O-P and P-O transitions 
are forbidden, has spectroscopic consequences. 



2.2 PES and the energy levels 

The energies of the quantised rotation-vibration states are 
the eigenvalue solutions that satisfy the Shrodinger equation 
for the oxygen and two hydrogen nuclei moving within the 
electronic potenial. However, the electronic potential within 
which the charged nuclei are moving is itself a function of 
the actual internuclear geometry. Therefore, in order to solve 
for the nuclear motion, it is necessary to have an accurate 
model of how the potential varies with the nuclear geom- 
etry. The problem is rendered tractable by adopting the 
Born-Oppenheimer approximation which separates the nu- 
clear and electronic motions and has as its basis the fact 
that, due to their lightness, the electrons may be considered 
to react immediately to any changes in the nuclear geom- 
etry. The most accurate PESs are computed using an ab 
initio starting point, with the resulting surface being em- 
pirically adjusted to improve the agreement between the 
computed energies and experimental data (Partridge and 
Schwenke 1997). 

We used the potential energy surface fit B of Shirin 
et al. (2003), which is based on a highly accurate ab ini- 
tio surface with adjustments for electronic relativistic and 
adiabatic (also known as the Born-Oppenheimer diagonal 
correction) effects and fitted to the available experimental 
data. At the time of writing, this surface is the most ac- 
curate available and it is the single most important factor 
affecting the accuracy of our results. 

Non-adiabatic corrections to the Born-Oppenheimer 
approximation are also important in the case of water 
(Schwenke 2003). However, a full theoretical non-adiabatic 
adjustement to the PES has been shown to produce no sig- 
nificant benefits compared to simplified approaches (Ten- 
nyson et al., 2002). These authors examine two such ap- 
proaches: using separate reduced masses for the vibrational 
and rotational motions (this systematically over-corrects 
bending motions and under-corrects stretches) , and a simpli- 
fied version of the full correction which includes only terms 
that scale with the kinetic energy terms in 8 and r. The latter 
approach is preferred and the ajustment is effected through 
the DVR3D program rather than changes to the PES. 
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2.3 Transition intensities 

The intensities of the allowed transitions between rotation- 
vibration states are determined by the dipole transition mo- 
ments for these pairs of states. They are: 

<V>W> (i) 

where tp' and ip" are the wavefunctions of the two states, and 
JI is the electronic dipole moment vector. Therefore, in order 
to calculate intensities, a dipole moment surface (DMS) is 
required. 

The parameter computed by the DVR3D program suite 
is the Einstein A coefficient for, Aj/, for each transition. This 
is the coefficient of spontaneous emission between the upper 
and lower states. It is related to the dipole transition mo- 
ment for the pair of states and to J for the upper state. Aj/, 
is independent of temperature, relates to a single molecule, 
has units of s _1 and is given by: 

A if = ^v 3 9iVJ" + mi>'\TiW')\ 2 (2) 

where prime and double prime represent the upper, i, and 
lower, f, states respectively. 

The quantity usually derived from observation, is the 
line intensity, I, which has units of cm/molecule. I is 
temperature-dependent and in emission is related to Aj/ 
by the expression: 

C(2J' + 1) ( ~hcE" \ r /-ftci/VI 
1 = Q.-flVg. 6XP {-kT- ) I 1 - ° XP {-W ) \ Af (3) 

where v is the frequency in cm -1 , E" is the energy of 
the lower ro- vibrational level in cm" . Q vr {T) is the ro- 
vibrational partition function and is dimensionless, and gi 
is the nuclear spin degeneracy and caries only one subscript 
since transitions between different nuclear spin states are 
not allowed (Miani and Tennyson 2004). Boltzmann's con- 
stant, k, has units of JK _1 and the constant C has the value 
(87TC)- 1 = 1.3271xl0" 12 s cm" 1 . 

Since the BT2 line list includes Aj/ for each transition, 
the above equation enables the line intensities to be com- 
puted at any given temperature. 

Unlike the PES, where the most accurate are ab initio 
surfaces that have been fitted to the available experimental 
data, in the case of DMS, limitations on the accuracy of ex- 
perimental line strengths means that at the current time the 
most accurate surfaces are purely ab initio (Lynas-Gray et 
al. 1995). Two such DMS were tested in the intensity part 
of our calculations. Our initial BT1 line list was computed 
using the same PES as BT2 and a preliminary Lynas-Gray 
ct al. (in preparation) DMS. The line positions in BT1 arc 
identical to those in BT2 (the same PES having being used 
for both). However, we observed that the Einstein A coef- 
ficients of the weaker lines generated using the Lynas-Gray 
et al. DMS were often too large when compared with ex- 
periment. This part of the computation was therefore re- 
peated using the DMS of Schwenke and Partridge (2000). 
The results, which are contained in BT2, show much better 
agreement between the computed strengths of weak lines 
and experiment. The Einstein A coefficients of the stronger 
lines in BT1 and BT2 generally agree to within two percent 
and both agree reasonably well with experimental values. In 
addition to being superior to the Lynas-Gray et al. surface 



that we tested, Schwenke and Partridge's DMS represents 
a major improvement on the earlier DMS of Partridge and 
Schwenke (1997) and from our analysis, as well from the 
same authors, is the most accurate in existence. 



3 CALCULATING THE LINE LIST 

The BT1 and BT2 line lists were computed using the 
DVR3D suite of programs (Tennyson et al. 2004) on three 
Sun 5 Microsystem V880 mainframe computers: Enigma and 
Ra which are clustered using high speed interconnects, each 
having 8 processors and 32 Gb of RAM and 432 Gb of disk 
storage and PSE, which has 24 processors and 96 Gb of 
RAM and 1,296 Gb of disk storage at UCL's Hiperspace 
computing centre. The final, 'DIPOLE' stage of the suite was 
amenable to parallelisation with little time penalty. Other 
parts of the program were run on single processors, which 
avoided coding problems and was more efficient in computer 
time. The total number of processor hours employed in gen- 
crating the BT1 line list (including the preliminary conver- 
gence testing which is discussed below) was 55,000 hrs and 
a further 10,000 hrs. were used in repeating the 'DIPOLE' 
runs to generate BT2 and in testing the outputs. 

DVR3D calculates the bound rotation- vibration energy 
levels, the wavefunctions on a grid in three dimensional 
space, and the dipole transition strengths of the allowed 
transitions. The final part of the DVR3D suite, 'SPECTRA' 
is able to compute temperature-dependent spectra over any 
selected frequency range, convolved with either the natu- 
ral line width or some other selected profile, such as that 
given by the resolving power of a particular spectrometer. 
DVR3D uses an exact (within the Born-Oppenheimer ap- 
proximation) kinetic energy operator. The program uses a 
discrete variable notation (DVR) with two radial and one 
angular co-ordinate for the nuclear motion problem (vibra- 
tional and rotational motions being treated spearately). 

In establishing the working parmeters for our calcula- 
tions, we aimed to provide a line list that would be complete 
and accurate at the temperatures of late series K-series stars 
(up to 4,000 K) and at wavelengths down to 0.8 fim. Pre- 
liminary calculations showed that in order to achieve this 
it would be necessary to include all states lying at energies 
up to 30,000 cm -1 relative to the ground state of the sys- 
tem. Previous workers (Miller et al. 1994 and Partridge and 
Schwenke 1997) selected lower energy cut-offs. 

Our cut-off for the total angular momentum was J=50. 
We calculate that the highest value of J which has ro- 
vibrational energies of less than 30,000 cm~l is 58. How- 
ever, we also estimate that by terminating our calculation 
at J=50 we omit less than 500 levels out of a total of more 
than 505 million and none of the omitted levels has an energy 
less than 23,490 cm -1 , with the majority being at energies 
above 28,000 cm -1 . Consequently, even at a temperature of 
4,000 K these missing levels contribute less than 0.02% of 
the total partition function. Moreover by omitting Js above 
50 we saved in the region of 8,000 processor hours. 

In order to generate accurate eigenvalues, it is essential 
that all calculations are fully converged within the limita- 
tions imposed by computing power and time. The lack of 
convergence in earlier lists has already been noted (Polyan- 
sky et al. 1997). 
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DVR3D has the option of using Jacobi or Radau co- 
ordinates. The latter was selected as being more appropri- 
ate for water. Our Radau grid is defined in terms of two 
radial co-ordinates, ri and r2, each with 28 points, and one 
angular co-ordinate, 9 having 44 points. These numbers are 
determined by convergence testing at J=20. The eigenvalue 
solutions are particularly sensitive to the number of radial 
points, r„. However, the computation time rises very rapidly 
as this number is increased, so the selection of r n represents 
a compromise and is the principal factor impacting on con- 
vergence. 

The angular grid points are arrived at using associated 
Legendre polynomials for the underlying basis sets, whilst 
the radial grid points are set up using Morse oscillator- 
like underlying basis sets, which are defined in terms of 
parameters, r e , uj e and D e . These must be entered into 
'DVR3DRJZ', the first module of the DVR suite. The pa- 
rameters have physical counterparts which are respectively: 
the equilibrium bond length, harmonic frequency and disso- 
ciation energy of the water molecule. However, the parame- 
ters in the Morse oscillator function differ from the physical 
values and must be determined by empirical testing. 

Although in principle, DVR3D is not strictly a varia- 
tional method, in practice it is found that the Variational 
Principle does apply; this fact is employed in obtaining val- 
ues for r e , L0 e and D e . The method is to alter the three 
parameters in a systematic manner until DVR3DRJZ gen- 
erates a set of pure vibrational states (the J=0, para states), 
the sum of whose energies is a minimum (251 eigenvalues 
were employed) . We conducted this part of the process man- 
ually. 

The investigation was complicated by the existence of 
local minima, which are not global minima for the three vari- 
ables. The set of parameters giving this result was: r e = 2.05, 
Lo e = 0.008 and D e = 0.20 all in atomic units (the equivalent 
dissociation energy is 43,895 cm -1 ). A good choice of these 
parameters (which are a function of the system and of the 
particular energy range that is of interest) is fundamental 
to the accuracy of the ensuing calculations. The eigenvalues 
are particluarly sensitive to r e and uj e and it was observed 
that differences of as little as 0.05 in the first of these two 
parameters and 0.002 in the second could affect the energy 
levels by in the region of 0.01cm -1 for states with energies in 
the region of 10,000 cm -1 . Larger deviations in these basis 
functions from their determined optimum values produced 
correspondingly greater errors and in the case of states with 
energies over 20,000 cm -1 , a bad choice of parameters could 
easily result errors in the individual levels in excess of 20 
cm -1 . Consequently considerable time was spent in deter- 
mining the values of the Morse oscillator- like basis set. 

Two other inputs are required by the vibrational mod- 
ule, DVR3DRJZ. These are: the maximum size of the inte- 
mediate Hamiltonian, which we chose after testing as 2500 
and the number of eigenvectors to be saved for use in the 
rotational module, 'ROTLEV3B', the optimum value for 
which was found by testing to be 700. 

ROTLEV3B also requies one variable to be determined; 
this is IBASS, the size of the Hamiltonian in the rotation 
module. IBASS varies with J. Its value was established by 
convergence testing at J=20 as being 530 x (J + 1 — p), 
where p is the rotational parity and has the value for even 
parity states and 1 for odd parity states. Hamiltonians of 
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A: Row in file, B: J, C: Symmetry (1-4: see Table 1), D: Row in 
block, E: v in cm -1 F, G, H: v\, U2, vo,. I, J, K: J, K a ,K c . 



this size are expensive in terms of computing time at high Js. 
Nevertheless, it was easy to demonstrate that lower values 
of IBASS produce results that are not converged. This is 
significant, as earlier workers using similar DVR techniques 
have used lower IBASS values. Viti et al. 1997, for example, 
used 200 x (J + 1 — p), and PS used an even lower effective 
number, particularly at high J. 

In selecting the various parameters referred to above, 
we regularly tested for convergence. We estimate that our 
choice of r n =28 accounted for approximately half of our total 
convergence error, which we estimated as being less than 
0.01 cm -1 at 10,000 cm -1 and in the region of 0.02 cm -1 
at 20,000 cm" 1 . 



4 RESULTS 

The BT2 water line list is available electronically in com- 
pressed form at: ftp://cdsarc.u-strasbg.fr/cats/VI/119 The 
data are in two parts. The first, the " Levels File" is a list 
of 221,097 energy levels, ordered by J and symmetry block. 
About 25,000 of these energy levels have been labelled with 
the appropriate angular momentum (3,K a ,K c ) and vibra- 
tional (^1,^2,^3) quantum numbers. An extract from the 
Levels File with an explanation of the contents of each of 
the 11 columns in the file is given in Table 2. 

The second part of BT2 is the "Transitions File". This 
has 505,806,202 entries. Each transition references upper 
and lower energy levels in the Levels File and gives the Ein- 
stein Aif coefficient for the transition. An extract from the 
Transitions File is given in Table 3. 

In uncompressed form the BT2 Transitions File is 12.6 
Gb of data. Therefore, in order to facilitate use of the list, the 
transitions have been ordered by frequency and separated 
into 16 smaller files, each representing a specific frequency 
range. 

In addition to the files containing the actual line list, 
the Strasbourg site contains a Fortran program, spectra- 
BT2.f90, that will enable users to generate emission or ab- 
sorption spectra from BT2 by specifying various parameters 
including: temperature, frequency range, cut-off intensity 
and line width. There is also a facility to generate spectra 
with full ro- vibrational assignments if required. The method 
of using spectra-BT2.f90 is detailed in a 'readme-spectra' file 
and there are also examples of a job file and an output file. 
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Table 3. Extract from BT2 Transitions File 
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A, B: Row numbers in the Levels File (upper and lower levels 
are not identified as the program tests for these). C: Aif (s _1 ). 



Table 4. Comparison of BT2 and PS (Partridge & Schwenkc 
1997) with 14,889 experimentally-determined energy levels 
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4.1 Comparing the BT2 and Partridge &: 
Schwenke line lists 

Several water line lists are in regular use by astronomers and 
the the most accurate list previously is that of Partridge and 
Schwenke (PS). Table 4 compares BT2 and PS energy levels 
with known experimental values (Tennyson et al. 2001). It 
will be seen from this table that although the PS list is more 
accurate than BT2 in the cases where agreement with exper- 
iment is better than 0.1 cm -1 , this is not the case generally. 
Specifically, based on a sample of 14,889 levels, whilst 99.9% 
of the BT lines are within 3.0 cm - of experiment 3.5% of 
the PS lines are outside this range. Other line lists (MT, 
VT2, SCAN) perform significantly worse than this. 

Examining deviations from experiment by energy is 
even more revealing, for it is seen that PS is increasingly 
unreliable above 10,000 cm - , which is the region with the 
greatest number of transitions. 

4.2 Labelling the levels 

As with observed water lines, the production of a synthetic 
line list prompts the question of how to assign quantum 
numbers to the transitions. The BT2 line list contains over 
505 million transitions, but the DVR3D suite only provides 
data on J and the symmetry block of each energy level. 
Since each individual line is a transition between two energy 
states, the problem reduces to one of labelling the 221,097 
energy levels, but this is still a large task. 

Many papers have been devoted to assigning quantum 
numbers to energy levels that have been deduced from ex- 
perimental line frequencies (e.g. Zobov et al. 2000). So far, 



Table 5. Distribution of levels in the BT2 and PS (Partridge & 
Schwenke 1997) disagceing with experiment by more than 2 cm -1 
- by frequency 
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1,914 








Total 


14,889 


30 


506 



less than 15,000 experimentally-determined levels have been 
labelled with their 3 rotational and 3 vibrational quantum 
numbers, even though much effort has been expended on 
the task. 25,870 levels in the BT2 list have been labelled 
using methodologies detailed below. BT2 labels 270 of the 
416 levels which have J=0 and energies below 30,000 cm -1 . 
However, the proportion labelled is less at higher Js. 

The process of labelling BT2 started by identifying par- 
ticular energies with experimental levels that have already 
been determined. In addition, we have labelled many of the 
energy levels that are unknown experimentally, using sev- 
eral different methods. A large number of vibrational states 
for J=0 have been labelled by visual inspection of the nodal 
structure of the wavefunction as described in Mussa and 
Tennyson (1998). This method has the disadvantage that 
observations can only easily be made on two-dimensional 
sections of the three-dimensional wavefunction and the pro- 
cedure can be misleading. In addition, Mussa and Tennyson 
oberved that at energies in the region of our 30,000 cm - , a 
high proportion of the wavefunctions are irregular with no 
identifyable nodal structure. . 

A second method involved examining the Aif coef- 
ficients for pure rotational transitions between an unla- 
beled vibrational state and a known vibrational state; the 
strongest transitions being those where standard selection 
rules are obeyed. We found this method (see Tolchenov et 
al 2005) to be useful for J<10 

A third method that was found to be useful in labelling 
higher J states involved the use of an algorithm to identify 
sets of levels within the same parity block having the same 
K a , z/i, V2, vz quantum numbers, but different values of J 
and K c . The method was originally developed for labelling 
the energy levels of the HCN and HNC isomers (Barber et al. 
2000). However, because of the density of the energy levels 
in water, an extra term was introduced when labelling the 
water levels. The algorithm used in this case was: 

E Jn £* AE Jn _ x - 6E Jn _ 2 + AE Jn _ 3 - Ej n _ 4 (4) 

However, when resonance between levels caused the be- 
haviour pattern to be erratic, or no Ej n _ 4 value existed, 
the original, simpler algorithm was used: 

E Jn £* ZEj^ - 3E Jn _ 2 + E Jn _ 3 (5) 

where Ej„ is the energy of the state in the same symme- 
try block having J = n and the same set of K , v\, 1/2, vz 
quantum numbers. 

The results of the labelling exercise are included in the 
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BT2 Level File. This means that when synthetic spectra are 
generated many of the transitions are fully labelled. This 
feature is useful when generating synthetic spectra for as- 
tronomical or laboratory applications as is discussed in the 
next section. 



4.3 Completeness 

If we compare the partition function, Q(T) for water com- 
puted at a particular temperature using the BT and PS line 
lists with the most accurately known value at this temper- 
ature, it is possible to estimate the the completeness of the 
line lists and the amounts of opacity that are missing in 
spectra generated by the two line lists at the selected tem- 
perature. 

A calculation of the partition function of water at 
3,000K using the 221,097 energy levels in BT2 yields a value 
that is 99.9915 percent of the Vidler & Tennyson (2000) 
value, which indicates that levels missing from BT2 only 
contribute about 85 parts in a million to the partition func- 
tion of water at this temperature, the reason being that there 
is a diminishingly small probability of states above 30,000 
cm -1 being occupied at this temperature. For comparison, 
the PS line list, which has 28,000 cm -1 cut-off gives a parti- 
tion function at 3,000 K that is only 99.493% of the Vidler 
& Tennyson value. 

Although the exclusion of levels above 30,000 cm -1 does 
not materially affect the completeness of the BT2 list, it does 
affect absorption at shorter wavelengths. If we consider a 
photon of of wavelength 1 pm (energy 10,000 cm - ). This 
is able to be absorbed by a water molecule in a particular 
rotation-vibration eigenstate provided that there is another 
eigenstate exactly 10,000 cm -1 above this lower state into 
which the molecule may be excited. It follows that since the 
BT2 line list has an upper cut-off of 30,000 cm -1 , none of 
the energy levels in the list above 20,000 cm" 1 are capable 
of being excited by a 1/itm photon, since there is no corre- 
sponding upper level. 

If we examine the extent to which Q(3,000K) computed 
from BT2, but excluding all levels above 20,000 cm" 1 , falls 
short of the Vidler & Tennyson value we will have an indica- 
tion (this is an upper limit as it takes no account of blending 
effects) of the opacity that has been excluded by adopting 
a 30,000 cm" 1 cut-off. Performing the calculation gives a 
shortfall of 0.83%. 

In the case of PS, only energy levels below 18,000 cm" 1 
are able to absorb a 1/im photon, and computing Q(3,000K) 
using only PS states up to 18,000 cm" 1 shows that these 
comprise only 98.37% of the Vidler value. Hence it will be 
seen that the opacity deficit at 1/xm is in the region of twice 
as great in the PS list as in BT and the ratio increases at 
shorter wavelengths. 



5 SAMPLE APPLICATIONS 

Although BT2 shows good agreement with experimentally 
known lines this is not a sufficient test of its accuracy as the 
PES used to generate BT2 was fitted to the known experi- 
mental data (this is also true of the PS line list). 

The most effective way of checking the accuracy of a line 



list is to test its ability to predict or identify previously un- 
known lines in astronomical or laboratory spectra. BT2 has 
been used successfully as outlined below. It should be noted 
here that some of the earlier spectroscopic applications used 
BT1 as they predate BT2. However, the line positions of the 
two lists are identical, and the intensities of the strong lines 
are similar for the two lists. 

5.1 Astronomical Spectra 

5. 1.1 e Indi Ba 

A synthetic spectrum generated at 1,500K using BT1 was 
able to reproduce the previously unknown absorption fea- 
ture observed at 1.554 /jm in the spectrum of the early T 
dwarf e Indi Ba, discovered by Volk et al. (2003), as being a 
blend of six water lines with no individual line contributing 
more than 25% of intensity (Smith et al. 2003). 

5.1.2 Comet 153P/Ikeya-Zhang (2002 CI) 

BT1 was used to compute the frequencies and Einstein A co- 
efficients of the 64 transitions (up to J=7) that make up each 
of the 7 hot bands of water detected in Comet 153P/Ikeya- 
Zhang. Dello Russo et al. (2004) applied these data in de- 
termining the rotational temperature of the comet on three 
dates. 

5.1.3 Temperatures of comet forming regions in early 
solar nebula 

The hot-band transitions identified by BT1 in Dello Russo et 
al. (2004) were classified into ortho and para using the sym- 
metry information contained within the BT2 energy file (see 
column C Table 1). Transitions between different nuclear 
spin parities can be ignored (Miani and Tennyson 2004). 
Dello Russo et al. (2005) were able to deduce the primordial 
O/P water composition of three comets: C/1999 HI, C/1999 
S4 and C/2001 A2 and hence the temperatures of the dif- 
ferent regions of the early solar nebula in which the comets 
were formed. The normal O/P ratio is 3:1, but since the 
lowest ortho level lies 23.8 cm - (~ 34 K) above the ground 
state (which is a para state), the O/P ratio is sub-normal at 
temperatures below ~ 50 K. A comparison of the Einstein 
A coefficients in BT2 with those actually used in BT1 shows 
that Dello Russo et al.'s results would have been the same 
had the later line list been used. 

5.1.4 Detection of water lines in nova-like object V838 
Mon 

Synthetic spectra generated by BT2 were used to identify 
absorption features observed in the 1.73 to 1.75 fim region 
of the spectrum of the nova-like object, V838 Monocerotis 
on five separate dates as being due to blended water lines. 
Quantum numbers were assigned to the 17 strong transi- 
tions that comprise the 5 absorption features (Banerjee et 
al. 2005). Sixteen of the lines were found to be in the (0 
0)-(0 1 1) band. This is the first time that individual water 
lines have been identified in a nova-like outflow region. 

In addition, BT2 was used to compute the theoretical 
intensities of the five absorption features as a function of 



BT2 water line list 7 



temperature and column density (assuming LTE). The re- 
sults indicated that the water features were arising from a 
cool ~ 750-900 K region around V838 Mon that was cooling 
at a rate of ~ 100 K per year. Column densities computed 
for the five dates also showed a reduction with time. 



5.1.5 Sunspot spectra 

Following their assignment of high temperature laboratory 
water lines in which extensive use was made of BT2 (see 
'Laboratory Spectra' below), Coheur et al. (2005) revisited 
the sunspot absorpton spectrum in the 9.89-12.95 /im region 
of Wallace et al. (1995) which had been partially assigned 
by Polyansky et al. 1997 and Zobov et al. 1999. Coheur 
et al. used their high temperature laboratory assignments 
to identify a substantial number of previously unassigned 
sunspot lines. 

5.2 Laboratory Spectra 

Hot water spectra have been analysed in the laboratory over 
the last thirty years. Methods generally involve high reso- 
lution Fourier transform spectrometry of vapour which may 
be at an elevated temperature, such as in an oxy-acetylene 
flame (Camy-Pyret et al. 1977, Coheur et al. 2005). 

The BT2 line list has already been used on a number 
of occasions to analyse spectra generated at both high and 
low temperatures at various wavelengths and this has has 
added considerably to the existing database of experimen- 
tally known energy levels and transitions. 

Coheur et al. (2005) working with experimental spec- 
tra generated at 3,000 K in the 5-20 /im region labelled 
about 600 previously unidentified levels using the BT2 line 
list. The identification of these levels was an important fac- 
tor in their subsequent assignment of 8,500 of the 10,100 
lines that they observed in the 0.385-1.081 /jm region. Most 
of the states that were labelled by Coheur et al. were ei- 
ther high J states having low bending modes or else lower 
J states with higher bending modes. These states had de- 
fied earlier analysis because the previous line lists used in 
earlier work, were unable to accurately predict the energies 
of states above 15,000 cm -1 or to treat accurately the high 
bending mode states, even those below 15,000 cm -1 . 

They note that whilst the agreement between BT2 and 
observation is generally within 0.1 cm -1 for the energy levels 
that they observe for low J states, the disageement can be 
as great as 0.8 cm -1 for some of the very high J states. Nev- 
ertheless, they comment that since for a given vibrational 
and K a state the diffence between the BT2 list and exper- 
iment increases smoothly with J, they are able to use BT2 
to predict the positions of unknown higher J levels with an 
accuracy of 0.02 cm -1 . This is considerably less than the 
experimentally determined line widths (full widths at half 
maximum) which were in the range 0.05-0.10 cm -1 due to 
broadening at 1 atm. and T = 3000 K. 

In a follow-up paper analysing the 2-5^im region Zobov 
et al., 2006 use BT2 to label approximately 700 previously 
unidentified ro- vibrational energy levels for water, observed 
in laboratory torch spectra. 

Tolchenov et al. 2005 analysed long-pathlength room 
temperature spectra. They use three separate line lists in 



their work. However, only the BT2 list is found to be reliable 
over all transition frequencies and in the case of lines having 
frequncies above 16,000 cm -1 it is the only one that they 
use. In addition, Tolchenov et al. encounter difficulties with 
the labelling adopted in previous studies, finding, for exam- 
ple, that in some cases different states have been labelled 
with the same quantum numbers. They therefore undertake 
a systemmatic re-labelling exercise using the BT2 energies. 
These labels are incorporated into the BT2 list for states 
with J < 9. 

Dupre et al.(2005) found BT2 similarly reliable for pre- 
dictions of long-pathlength room temperature spectra in the 
near ultra-violet. They observed 62 R-brach transitions in 
the 8v polyad and were able to determine 36 energy levels, 
previously unknown experimentally. 



6 CONCLUSION 

We present a new synthetic water line list which gives the en- 
ergies of 221,097 states with cut-offs of J=50 and E=30,000 
cm -1 . 25,870 of the lower energy levels have been labelled 
with a full set of three rotational and three vibrational 
quantum numbers. BT2 lists 505,806,202 trasitions. It has 
been extensively tested against experimental observations 
and also compared with other lists. It has been shown to be 
the most complete and accurate water line list in existence. 

We make our results freely available in electronic form 
via ftp://cdsarc.u-strasbg.fr/cats/VI/119, in the hope that 
BT2 line list will be a vauable tool for astronomers in both 
spectroscopy and atmospheric modelling applications. 

One of the problems facing the modellers of the atmo- 
spheres of cool stars and brown dwarfs is the disageement 
between observation and model. Nevertheless, considerable 
progress has been made in recent years in such areas as: 
convection, molecular abundances, non-LTE effects and dust 
(Hauschildt et al. 1999; Tsuji 2002; Alexander et al. 2003). 
At the same time, there have been advances in computing 
the opacity effects of the various species that are included 
in the models. 

At the temperature of late M-dwarfs (2,500K), water is 
by far the most important contributor to stellar atmospheric 
opacity, typically contributing over 60% of all opacity in the 
infra-red. It is also an extremely important contributor to 
the opacity of L and T brown dwarfs, although at these lower 
temperatures (down to 900K), methane and ammonia play 
an increased role. 

As part of our on-going work in applying BT2 we plan 
to compare observed M-dwarf spectra with synthetic spec- 
tra produced with the Phoenix model (Hauschildt, Allard 
& Baron 1999) using several water line lists, including BT2. 
This is an extension of work already conducted on mod- 
elling oxygen rich cool stars. Jones et al. (2005) compare 
the CO 2-0 bands in the 2.297-2.210 ^m region of several M 
dwarfs and an L dwarf with opacity calculations using the 
PS and BT line lists and conclude that for this particular 
wavelength range 'While the Partridge- Schwenke line list is 
a reasonable spectroscopic match [for the BT line list] at 
2,000K, by 4,000K it is missing around 25 % of the water 
vapour opacity.' 
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